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ABSTRACT 

In this paper, we have presented a method of evaluation of critical 
frequent Q c of vortex whose circulation number k=1 confined in a 
spherical trap with N atoms of 87 Rb. We have taken the values of 
a h0 =0.791 um. Our K=1 theoretically evaluated results indicate that 
Q c decreases with N and decrease is more pronounced as one goes 
to higher N. Our theoretical results are consistent with those of the 
others theoretical workers. 
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1. INTRODUCTION 

Super fluidity is one of the most 
spectacular consequences of Bose-Einstein 
condensation. However, the explicit 
connection between super fluidity and BEC 
is not trivial and has been the object of a 
long-standing and deep investigation in the 
last decades, mainly forits importance in 
understanding the physics of liquid helium. 
In macroscopic bodies super fluidity shows 
up with many peculiar features : absence of 
viscosity, reduction of the moment of inertia, 



occurrence of persistent currents, new 
collective phenomena (second sound, third 
sound, etc.), quantized vortices and others. 
Several properties are usually interpreted as 
coherence effects associated with the phases 
s of the super fluid through u s =(h/m) V s. A 
major question is to understand whether 
some of these effects can be observed also in 
trapped gases. Of course in a mesoscopic 
system, one expects the manifestations of 
super fluidity to be different from the ones 
exhibited by macroscopic bodies. In 
particular, traditional experiments based on 
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the study of transport phenomena are not 
easily feasible in trapped gases. On the other 
hand, interference patterns, associated with 
phase coherence, have been already 
observed (Andrews, Townsend et.al.) 1 and 
successfully compared with theory. This 
opens a promising field of research based on 
the investigation of coherence phenomena, 
including the realization of the so-called 
"atom laser." 

2. MATHEMATICAL FORMULA USED 
IN THE EVALUATION 

Among the several properties 
exhibited by super fluids, the occurrence of 
quantized vortices and the strong reduction 
of the moment of inertia represent effects of 
primary importance. 

In a dilute Bose gas, the structure of 
quantized vortices can be investigated 
starting from the Gross-Pitaevskii equation. 
Indeed one of the primary motivations of the 
theory was the study of vortex stated in 
weakly interacting bosons. 2 4 These studies 
were further developed by Fetter, 5 including 
higher-order effects in the interaction. 

A quantized vortex along the z axis 
can be described by writing the order 
parameter in the form 



0(r) = 0 v (r ± ,z)exp(ixp) 



(1) 



where (p is the angle around the z axis, K is 

an integer, and 0 V (r ± , z) = yj n(r ± , z) This 
vortex state has tangential velocity 



v = - 



K 



mr. 



(2) 



The number A" is the quantum of circulation 
and the angular momentum along z is NKh . 
The equation for the modulus of the order 
parameter is obtained from the GP equation. 
The Kinetic energy brings a new centrifugal 
term arising from the velocity flow that 
pushes the atoms away from the z axis. The 
GP equation then takes the form 



r n 2 v 2 s¥ m, , 2 , . 2 , 



Ki,z) = 



Ul, z) 
(3) 



Due to the presence of the 
centrifugal term, the solution of this 
equation for K has to vanish on the z axis. 
For non interacting systems the solution of 
Eq. (3) is analytic and for K-\, has the form 



^(r J _,z)«,r J _exp 



mi 2 2 \ 



(4) 



In this case, the vortex state 
corresponds to putting all the atoms in the 
m=l single-partial state. Its energy is then 
NhC0 ± plus the ground state energy. 

Similarly to what happens for the ground 
state without vortices, the presence of 
repulsive forces; dramatically reduces the 
density with respect to the non interacting 
gas, the condensate wave function becoming 
much broader. 

The structure of the core of the 
vortex is fixed by the balance between the 
kinetic energy and the two-body interaction 
term. For a uniform Bose gas, the size of 
the core is of the order of the healing length 

E,={%7lna) " 2 , where n is the density of 

the system. For the trapped gas, an estimate 
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of the core size can be obtained using the 
central value of the density in the absence of 
vortices for n. If the trap is spherical, the 
ratio between and the radius R of the 
condensate takes the form 6 . 



I 
R 



( . 



"ho 

R 



(5) 



where we have used the Thomas-Fermi 
approximation for the central density and the 
radius R. For the condensate in the figure, 
the radius is about 4.1 in units of a ho and the 
ratio £ I R is then ~ 0.06. The actual core 
size depends, obviously, also on the position 
on the position on the z axis and becomes 
larger when the vortex line reaches the outer 
part of the condensate, where the density 
decreases. 

The energy of the vortex can be 
evaluated through the energy functional (5). 
The difference between the energy of the 
vortex state and the one of the ground state 
allows one to calculate the critical frequency 
needed to create a vortex. In fact in a frame 
rotating with angular frequency Q. the 
energy of a system carrying angular 

momentum L z is given by (E-£lL z ), 

where E and L z are defined in the laboratory 
frame. At low rotational frequencies this 
energy is minimal without the vortex. If D. 
is large enough the creation of a vortex can 
become favorable due to the term -Q.L, . 
This happens at the critical frequency 



fl, = (^[(-),-yj 



(6) 



where E K is the energy of the system in the 
presence of a vortex with angular 



momentum NfiK . When one plots the 
critical frequency for the creation of a vortex 
with K — 1 as a function of the number of 
atoms, for rubidium in a spherical trap. As 
shown in the figure, the predicted critical 
frequency is a fraction of the oscillator 
frequency and, in typical experimental 
conditions, corresponds to a few Hz. It 
decreases when N increases, because for 
large systems the energy cost associated 
with the occurrence of vortex increases as In 
N, while the gain QL z is always linear in N. 

This behavior is similar to the one exhibited 
by uniform systems where, approximating 
the vortex core as a cylindrical hole of radius 
i^, one can explicitly calculate the critical 

frequency and finds Q. l = (N I mR 2 ) in 

[R I , where R is the radius of the region 

occupied by the vortex flow. Analogous 
expressions can also be derived in the 
presence of harmonic trapping for large N. 
Baym and Pe thick 6 and Sinha 7 have shown 
that the critical frequency, in units of CO ho 

goesD (a ho I R) 2 in{RI^). Using the 

asymptotic solution of the GP equation in 
the large N limit, Lundh, Pethick and Smith 8 
have found a useful analytic expression for 
the critical velocity in the case of axially 
symmetric traps: 



£1 = 



5h , 0.6717? 

-In - 



2mR\ 



(7) 



Where R ± is the Thomas-Fermi radius of 

the cloud in the xy plane orthogonal to the 
vortex line, while the healing length is 

defined by B,-{%7tna) '"with n equal to 
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the central density of the gas without vortex. 
This formula gives a critical frequency that 
significantly differs from the numerical 
result N smaller that about 2000, while for 
larger N it becomes more and more accurate. 

3. DISCUSSION OF RESULTS 

In this paper, we have presented a 
method of evaluation of critical frequency 
Q. c of vortex whose circulation number 

K = 1 confined in a spherical trap with N 
atoms of we have taken the value of 
a ho =0.791 jUm.In table Ti we have given the 

evaluated results of critical frequency Q. c in 

the units of Q) L as a function of atoms N. We 

have also compared our theoretically 
evaluated results with those of the other 
worker. 9 14 

Table T t 

Evaluated results of critical angular velocity 
Q. c in units of CO L for the vortex K — 1 in a 
spherical trap as a function of atoms N of 
87 s The value of a h =0.791//m 





Q. c in units of 0) ± 


N 


Ours 


Other 6 


0 


1.05 


1.16 


1000 


0.85 


1.10 


2000 


0.82 


0.96 


3000 


0.78 


0.82 


4000 


0.68 


0.76 


5000 


0.46 


0.70 


6000 


0.42 


0.62 


7000 


0.40 


0.48 


8000 


0.38 


0.42 


9000 


0.36 


0.40 


10000 


0.35 


0.40 



Evaluated results of (© / & ri ) as a function 

of (TIT 0 ). 

(© / ® rig ) = Interacting gas with 7] =0.4 
(0 / 0 nt ) 2 = non-interacting gas with 

N=5Xl0 7 atoms 
(0/0 r)g ) = non-interacting gas with 

N=5Xl0 4 atoms 



(t/t c °) 


(© /©,,)' 


(0/©„ s ) 2 


{e/e ris f 


0.1 


0.052 


0.087 


0.142 


0.2 


0.087 


0.112 


0.262 


0.3 


0.104 


0.147 


0.305 


0.4 


0.118 


0.236 


0.386 


0.5 


0.287 


0.305 


0.408 


0.6 


0.326 


0.368 


0.435 


L 0.7 


0.453 


0.488 


0.552 


0.8 


0.567 


0.605 


0.786 


0.9 


0.669 


0.708 


0.835 


1.0 


0.728 


0.814 


0.926 



Q. c Decreases with N and decrease 

is more pronounced as one goes to higher N. 
Our theoretical results are consistent with 
those of the other 6 as far as the trend is 
concerned. In other calculation, we have 

evaluated the ratio of (0/0^) as a 

function of {T IT®}. These results are 

shown in table T 2 .Here we have evaluated 
the ratio in three different cases. In first case 
we have taken the interacting gas with 
J] = 0.4 and in other two cases we have 
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taken the non-interacting gas having no. of 
atoms 5xl0 7 ara/5xl0 4 and respectively. 
In all the three evaluated results (@/0 n? ) 

increases with (T I T° ). 

The increases in more faster in the 
non-interacting gas case and less faster in 
interacting case. 15 Some recent calculations 
16 23 also revel the same behavior. 
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